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(N ; ABSTRACT 

1^ ' We present evidence for shock acceleration of cosmic rays to high energies (~10 TeV) 

G^ ■ in the supernova remnant IC443. X-ray imaging spectroscopy with ASCA reveals two 

regions of particularly hard emission: an unresolved source embedded in an extended 

Ph| emission region, and a ridge of emission coincident with the southeastern rim. Both 

P ■ features are located on part of the radio shell where the shock wave is interacting with 

molecular gas, and together they account for a majority of the emission at 7 keV. Though 

c^ ! we would not have noticed it a priori, the unresolved feature is coincident with one re- 

solved by the ROSAT HRI. Because this feature overlaps a unique region of flat radio 
spectral index [a < 0.24), has about equal light-crossing and synchrotron loss times, and 
a power law spectrum with a spectral index of a=1.3±0.2, we conclude that the hard 

p^ ', X-ray feature is synchrotron radiation from a site of enhanced particle acceleration. Evi- 

dence against a plerion includes a lack of observed periodicity (the pulsed fraction upper 
limit is 33%), the spectral similarity with the more extended hard region, the location 
of the source outside the 95% error circle of the nearby EGRET source, the fact that it 
is nestled in a bend in the molecular cloud ring with which IC443 is interacting, and the 
requirement of an extremely high transverse velocity (>5,000 km/s). We conclude that 
the anomalous feature is most likely tracing enhanced particle acceleration by shocks 
that are formed as the supernova blast wave impacts the ring of molecular clouds. 

Subject headings: acceleration of particles — radiation mechanisms: non-thermal — 
shock waves — cosmic rays — supernova remnants — supernovae: individual (IC443) 
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1. Introduction 

It is generally accepted from energy budget con- 
siderations that supernova remnants (SNRs) are the 
primary source of galactic cosmic rays (CRs), which 
are inferred to span the energy range from about 10^ 
to 10" eV (Blandford & Ostriker 1978; Axford 1994; 
Biermann 1995, and references therein). Though 
radio observations of synchrotron radiation from su- 
pernova remnants supply bountiful evidence that they 
are the source of CRs at GeV energies, evidence re- 
garding higher energy cosmic rays is largely circum- 
stantial. A search for ultra-high energy 7-rays ('^lO^'* 
eV) coincident with SNRs detected nothing signifi- 
cant (Allen, et al. 1995). The best evidence thus far 
for shock acceleration of ^100 TeV CRs, comes from 
X-ray observations of SN 1006 (Koyama et al. 1995). 
SN 1006 appears to be unique in that the non-thermal 
X-ray flux from its shell dominates the thermal, but 
current limits on the non-thermal flux from the shells 
of other young SNR do not exclude the presence of 
a synchrotron component, and hence ongoing CR ac- 
celeration to TeV energies. 

In this paper we report the discovery of a site in 
another SNR in which high energy cosmic rays are 
possibly being accelerated, but by a different mecha- 
nism from that in SN 1006. Using ASCA CIS data, 
we have found a concentration of hard X-ray emission 
along the southern rim of the middle-aged remnant 
IC443, whose spectrum is consistent with a power 
law, and whose morphology suggests enhanced shock 
acceleration resulting from the SNR shock encoun- 
tering dense clouds in the ISM. The observations are 
consistent with the prediction of such an effect by 
Jones & Kang (1993, hereafter JK93). 

IC443 (G189.1-f3.0) is a nearby (« 1.5 kpc) super- 
nova remnant (SNR) of intermediate age. Many in- 
frared and radio line observations have demonstrated 
that IC443 is interacting with a ring of molecular 
clouds (e.g., DeNoyer 1979; Burton et al. 1988, here- 
after BGBW; Burton et al. 1990; Dickman et al. 1992, 
hereafter DSZH; van Dishoeck et al. 1993). The re- 
gion in which the most complex interactions between 
the molecular cloud and the SNR shock is occurring 
also has an unusually flat radio spectral index for a 
shell-like SNR (a <0.24 - Green 1986, hereafter G86). 

In the X-ray band, IC 443 has been the subject of a 
number of comprehensive studies, most notably those 
of Petre et al. (1988, hereafter PSSW) and Asaoka & 
Aschenbach (1994, hereafter AA94). It has an irregu- 



lar soft X-ray morphology, influenced strongly by its 
interactions with an H I cloud to the northeast and 
the foreground molecular cloud and shadowing by the 
foreground SNR G189.6-h3.3. Using Gmga, Wang et 
al. (1992, hereafter WAHK) resolved spectrally a hard 
component which they were able to characterize by ei- 
ther a thermal model with kT~14 keV or a power law 
with photon index T ~2.2. Scanning observations by 
HEAOl A-2 showed that the hard emission is more 
centrally located than the bright soft emission. Nei- 
ther the A-2 nor the Ginga observations provided in- 
formation regarding the extent of the hard emission. 
IC 443 is also coincident with a high energy 7-ray 
source, which has led to speculation that the 7-rays 
are produced by the interaction of cosmic rays, ac- 
celerated by IC443's shocks, with nearby molecular 
cloud material (Sturner & Dermer 1995; Esposito et 
al. 1996, hereafter EHKS). Models of broad-band 
non-thermal spectra have recently been produced by 
Sturner et al. (1996); in these synchrotron radiation 
dominates from the radio to the soft X-ray, while elec- 
tron bremsstrahlung and tt'^ decay dominate the 7-ray 
emission for a supernova remnant like IC443. 

Our paper is organized as follows. We first present 
the ASCA GIS images showi ng the location of the 
hard X-ray emission region 



§2.1 



We compare the 
ASCA image with that from the ROSAT HRI. We 
then discuss our analysis of the GIS spectrum and 



the GIS and HRI light curves from that region (§2.4). 
Finally, we discuss possible emission mechanisms and 
the implications for the acceleration of cosmic rays by 
supernova remnants (§ [3.2| ). 

2. Observations and Analysis 



Table 1: ASCA observations of the SNR IC443 



Sequence # 


10011010 


51023000 


Mission Phase 


PV phase 


AO-1 


RA (J2000) 


06'^17'"25" 


Qgh^gm^gs 


Dechnation 


+22°40'37" 


+22°26'53" 


GIS Exposure 


20.2 ks 


34.7 ks 


Date 


April 14, 1993 


March 9, 1994 



2.1. The ASCA GIS X-ray Maps 

ASCA performed observations of two adjacent re- 
gions of IC 443, one during the Performance Veriflca- 
tion (PV) phase of of the mission, the other during the 
first cycle of guest observations (AO-1). The relevant 




Fig. 1. — Soft (A) and hard (B) GIS contour images 
of IC443. The energy bands are 1.1-2.1 and 2.1- 
12.0 keV respectively. The ratio of the two bands (C) 
is also shown (blanked in regions of low soft X-ray 
surface brightness). The contour levels in each figure 
are multiples of 10% of the peak. For reference, we 
have shown in grey-scale (A, B, C), a 1.4 GHz radio 
map from the VLA sky survey (10-100 mJy/beam). 
Image D shows our ROSATRRl mosaic of IC 443; the 
GIS field of view for each observation and the region 
shown in Figure are overlaid. 



information about these observations is contained in 
Table m We extracted these data sets from the 
ASCA public archive. In figures [^ A and B, we show 
exposure-corrected GIS mosaic images for the hard 
(2.1-12 keV) and soft bands (1.1-2.1 keV). These 
appear highly correlated, except for a bright feature 
in the hard band map, centered at RA ~ 06''17™05'', 
Dec = -h22°21'30" (J2000). The anomalous nature 
of this feature is shown in dramatic fashion in figure 
yC, a spectral hardness ratio map, constructed by di- 
viding the hard band map by the soft. There it can 
be seen that the other features, associated with the 
diffuse thermal emission in IC 443, all have similar 
spectral hardness, but the hard X-ray feature (HXF) 
stands out. The ASCA GIS brightness distribution 
within the hard feature is consistent with that ex- 
pected from a source smaller than ^ 1'. The hard 
feature is ~12' off axis so it was not in the SIS field 
of view. 
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Fig. 2. — The ASCA GIS hardness ratio map 
{F^.i-ii kev/^i. 1-2.1 kcv) is shown in grey-scale. The 
overlaid drawings are: the EGRET 7-ray detection 
circle (EHKS), the region of flat radio spectral index 
(G86) and a contour plot of the v=l-OS(l) II2 emis- 
sion line intensity (BGBW). 



In addition, there is a ridge of spectral hardness 
located northeast of the HXF along the radio-bright 
rim of the remnant. This second feature is not shown 
in Fig. I^C because it was blanked due to its small 
soft X-ray count rate. It is shown, however, in Figs. 
H and ||. As will be discussed in more detail in ijpii , 



the hard flux from the observed portion of the ridge 
is about half the hard flux from the HXF and the two 
regions can account for a majority of the total hard 
ASCA flux. 

2.2. i?05.4THRI Image 

The region containing the HXF was observed for 30 
ks using the ROSATYiYd as part of a program to cre- 
ate a complete high resolution X-ray image of IC 443. 
In Fig. np, we show the complete HRI mosaic, with 
the GIS fields of view overlaid. A surface brightness 
enhancement appears approximately at the location 
of the HXF. The HRI map of the region containing 
the HXF is shown in Fig. H, with the GIS hardness ra- 
tio contours overlaid. The feature's central core has 
an angular extent of about 10" x 5", with extended 
lower level structure of about 1' in extent. The elon- 
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Fig. 3. — The regions where ASCA GIS data were 
selected for our spectral analysis and the approximate 
ASCA GIS fields of view, overlaid on the same hard- 
ness ratio map as shown in Fig. ||. We purposely chose 
a larger than needed region about the HXF and the 
ridge in order to spectrally distinguish between the 
thermal and non-thermal spectra. 



gation of the central core is perpendicular to the ori- 
entation of the low level emission and does not resem- 
ble the HRI point spread function, so we believe it to 
be resolved. 

2.3. Spatial Correlations with Other Bands 

In Fig. g we overlay on the GIS spectral hard- 
ness ratio map possibly related spatial information 
from other bands: the EGRET 95 percent confi- 
dence error circle (EKHS), contours of II2 emission 
(BGBW) which locate the sites of the most intense 
SNR shock/cloud interaction, and the region in which 
the AA20-200 cm spectral index is appreciably flatter 
than elsewhere in IC 443 (G86). The HXF is just out- 
side a concave arc of H2 emission, but correlates well 
with the region of flat radio index. Interestingly, it 
is well outside the EGRET error circle. In the Digi- 
tized Sky Survey (Lasker et al. 1990), there are some 
faint filaments in this region, but no obvious unre- 
solved sources near the hard feature. A search of the 
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Fig. 4.— The 10" FWHM Gaussian smoothed 
i^OS'^T HRI (grey-scale) and the 60" resolution hard- 
ness ASCA hardness ratio map shown with the same 
contour levels as figure |l|C. The grey-scale range is 
10-45 counts per 10" diameter circular "beam." 



most recent on-line pulsar catalog (Taylor et al. 1993) 
reveals none near the feature. 

In addition to the main hard feature, the ridge of 
hard emission overlaps the H2 emission region, but 
does not appear coincident with any region of radio 
spectral index flattening. However, it does lie along 
the bright radio rim. 

2.4. The ASCA Spectrum 

The HXF and the ridge are embedded within the 
diffuse emission of IC443 — at least in projection. 
Thus in order to determine the spectral properties of 
these features, one must first adequately character- 
ize the "background" thermal emission. We there- 
fore extracted a spectrum from the adjacent larger 
region shown in figure 0. We assume, based on the 
smoothness of the spectral hardness ratio map, that 
the spectrum of the diffuse emission does not vary 
appreciably in the neighborhood of the HXF; based 
on the previous X-ray observations the spectral vari- 
ations in the diffuse emission would be most signifi- 
cant below 1 keV, outside the effective GIS band. To 
this "background" spectrum we fit an ad hoc thermal 
emission model, comprised of a thermal bremsstrah- 
lung continuum and narrow Gaussians to represent 
the most prominent emission lines (the Hea transi- 
tions of Ne, Mg, Si, S and Ar). Such ad hoc models 
have been used previously to model ASCA spectra 



Table 2: Spectral Fits to Backg 


round Field and Hard X-ray 


Features 








Model 


Data Sets 


x' 


xi 


Tsolt 

keV 


1022cm-2 


Hard 
Component 


1 


Thermal" 


Background 


781 


0.71 


0.89 ± 0.03 


0.12 ± 0.03 


n/a 


2 


Thermal" 


OnHXF 


793 


0.73 


1.9 ±0.1 





n/a 


3 


Thermal^ 


OnHXF 


1407 


1.28 


0.89 


0.31 ± 0.03 


n/a 


4 


Power-law'^ + Thermal'' 


On HXF 


769 


0.70 


0.89 


0.18 ± 0.03 


r = 2.3 ±0.2 


5 


Power-law + Thermal'^ 


On HXF 


766 


0.70 


-, 1+0.3 
-■^■^-0.1 


0.12 


r = 2.4io-i 


6 


Hot Bremss + Thermal'' 


On HXF 


763 


0.70 


0.89 


0.13 ± 0.03 


kT=3.9J:J-6 keV 


7 


Thermal"* 


The Ridge 


778 


0.40 


0.74 ± 0.08 


0.12 


n/a 


8 


Power-law -I- Thermal'' 


The Ridge 


755 


0.39 


0.89 


Q<0.03 


r = 0.1 ±1.3 


9 


Power-law -I- Thermal"* 


The Ridge 


748 


0.38 


^4+"-07 

^■3^-0.16 


0.12 


r = i.Qtl', 



Fits are to the energy range 0.5-12 keV 

" Letting the scaling (GIS2 & GIS3), column density, temperature and line strengths vary 

** We froze the line strengths and temperature to the off source (model #1) values, but fit the column density and 

scaling. 

"" The best-fit 5 keV fiux density is (7 ± 3) x 10-i3_£^, 

'^ We froze the column density and line strengths to the off source (model #1) values, but fit the temperature and 

scaling. 
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Fig. 5.— The ASCA GIS2 spectrum of the hard 
feature in IC443 with the best-fit ad hoc thermal 
model. By adding a hard power law component to 
this model we reduce the combined GIS2 and GIS3 
X^ by Ax^ = 640. 



from other remnants (e.g. Miyata et al. 1994; Holt 
et al. 1994). The best fit model yields kT~0.9 keV, 
comparable with published temperatures for IC 443 
(PSSW, WAHK). However, it is important not to 
over-interpret our thermal results, because we have 
made no attempt to distinguish the Hea transitions 
from weaker nearby spectral lines (including each el- 
ement's corresponding H-like Lya transition) or to 
develop a physically self-consistent thermal model. 



We next applied this same model to the GIS spec- 
trum of the region containing the HXF (also shown 
in figure 0), allowing only the normalization to vary. 
As shown in table p| this yields an unacceptable fit. 
As indicated in FigT]5[ the "background" model pro- 
vides a good fit to the data below ~2 keV, but there 
is a substantial flux excess at higher energies. An ac- 
ceptable fit to the "on HXF" spectrum is obtained 
by adding a second continuum component. The cur- 
rent data do not allow us to distinguish between a 
power law with r=2.3±0.2 and thermal bremsstrah- 
lung with kT=3.9lo:6 keV. 

We also performed this same analysis using a 
smaller on-source region, tightly drawn around the 
HXF. However, this was less effective: because of 
the lack of thermal background counts, it is harder to 
spectrally distinguish the background from the non- 
thermal source emission. Nevertheless we measure a 
photon index of r=2.2^Q'4, using the technique de- 
scribed above. 

We applied the same spectral analysis procedure 
as described above for the HXF to the hard ridge. 
Since the ridge is farther from the "background" re- 
gion than the HXF (see Fig. y) and we fit the ther- 
mal model using the AO-1 data, we were uncertain if 
it was appropriate to use the same background ther- 
mal model for both the HXF and the ridge. To test 



this, we performed our analysis using both the AO-1 
"background" region shown in Fig. and one using 
PV phase data closer to the ridge. We do not detect 
any significant difference in the spectral index or flux 
of the hard ridge region as a function of the chosen 
background region. 

Because of its larger extent and lower surface 
brightness, the ridge spectrum has lower signal to 
noise than the HXF. As with the HXF, fitting the 
spectrum with the background model yielded a clear 
excess above 2 keV. It was not possible, however, to 
constrain the spectral properties of the ridge's hard 
component. While we have no evidence suggesting a 
different spectral form, we cannot eliminate the pos- 
sibility that the hard emission from the ridge has a 
different origin. 

The GIS spectral parameters for this region are 
similar to those obtained using the non-imaging (l°field 
of view) Ginga LAC. WAHK reported the presence 
of a hard component that was fit by a thermal model 
with kT^ll^jQ kcV (from their Fig. 3) or a power- 
law with photon index F = 2.2 ± 0.13. Assuming a 
uniform surface brightness inside the shell and an ab- 
sorbing column density of A^h = 10^^'^ cm^^, WAHK 
found the total 2-20 keV emitted flux from IC443 
to be 9 X 10~^^ ergcm"^ s~^. As the total flux is 
dominated by the ~ 1 keV thermal component, a 
comparison between the broad band Ginga flux and 
the GIS flux from the hard features is inappropri- 
ate. A more reasonable quantity to compare is the 
flux density at 7 keV, where the relative contribu- 
tion of the ~ 1 keV component is negligible. Using 
Figure 2 in WAHK, and assuming that the Ginga 
LAC has an effective area of 3000 cm^ and near 
unity quantum efficiency at 7 keV, we estimate the 
Ginga flux density at 7 keV to be 20 x 10~^ pho- 
tons s^^ keV^^ cm^^. From the GIS data, we find 
the HXF flux density at 7 keV to be approximately 4 
X lO^'"^ photons s~^ keV~^ cm~^, and the flux den- 
sity of the portion of the hard ridge within the field of 
view to be about 2 x 10~^ photons s~^ keV""'^ cm~^. 
Thus the detected flux in the hard regions seems to 
account for only 30% of the hard flux detected by 
Ginga, despite their prominence in the GIS images. 
Moreover, integration over the rest of the surface of 
the remnant contained in the GIS fields (about 90 
percent of the total) yields a fiux density of no more 
than 4 x 10~^ photons s~^ keV~^ cm^^, leaving a 
factor of two discrepancy between the GIS and Ginga 
7 keV flux densities. 



Two possible resolutions of this discrepancy are 
that there is signiflcant hard flux from the small frac- 
tion of IC 443 unobserved by the GIS, or there is a 
calibration discrepancy between the GIS and Ginga. 
To check this latter possibility, we compare the broad 
band fiux from all of IC 443 in the common 2-10 
keV band. The 2-20 keV Ginga flux is 9 x 10"" 
ergs cm-2 s'^ (WAHK). The 2-10 keV fiux from IC 
443 as measured by HEAOl A-2 is 7±1 x 10"" 
ergs cm~^ s"^ (PSSW), consistent with the Ginga 
measurement. Integrating the 2-10 keV flux within 
the field of view of the two GIS pointings, we find 5±1 
X 10"^^ ergs cm"^ s"^; this represents 90 percent or 
more of the total fiux from the remnant. The con- 
sistency of these numbers indicates that contribution 
of cross calibration uncertainties to the discrepancy is 
probably small, but could be as large as a factor of 2. 
Thus while it is possible that the GIS has observed all 
the regions in IC 443 producing hard fiux, we cannot 
rule out the possibility that the unobserved portion 
of the IC 443 rim to the south and east of the ridge 
contribute significantly to the hard flux. 

A key conclusion about the spatial distribution of 
the hard X-ray emission in IC 443, independent of a 
resolution of this discrepancy, is that the hard emis- 
sion is highly localized. The fiux density from the bulk 
of the remnant surface area is at most 40 percent of 
the total hard flux. This provides a strict upper limit 
on the luminosity of a hard thermal component, dis- 
tributed throughout the remnant. WAHK inferred 
that IC443 is a very young (1200 y) remnant based 
on their interpretation of the hard component as hot 
gas arising from a high velocity shock, and IC443's 
coincidence with the guest star of 837 A.D. Despite 
the restrictive limit we have placed on the fiux of a 
hot thermal component, we nevertheless cannot rule 
out their conclusion about the age of IC443. More- 
over, high shock velocities may also be required to ac- 
celerate electrons to the relativistic energies required 
for the emission of X-ray synchrotron radiation, and 
WAHK's historical evidence is naturally unaffected 
by our observations. 

2.5. Pulsation Search 

In order to search for periodicity from the un- 
resolved HXF, we performed coherent fast Fourier 
transforms (FFTs) on the ROSAT HRI and ASGA 
events recorded from the region used for spectral anal- 
ysis (after applying the standard barycenter correc- 
tions). All our results are consistent with a non- 



periodic signal. 

The temporal resolution of the GIS is 0.0625 s, so 
we cannot use the GIS data to search for a frequency 
faster than 8 Hz. On the other hand, the ROSAT 
HRI data have millisecond resolution. We therefore 
performed summed FFTs on the combined data set 
{ROSAT + ASCA) in the low frequency range 0.05-8 
Hz, and only the ROSAT data, in the high frequency 
range 10-1000 Hz. We find 99% confidence pulsed 
fraction upper limits of ^15% and '^33% in the low 
and high frequency ranges respectively. 

3. Discussion 

3.1. Summary of Results 

We have discovered an isolated hard X-ray emit- 
ting feature and a ridge of hard emission in the south- 
east of the SNR IC 443 using the ASCA GIS. The 
HXF's X-ray spectrum can be characterized by ei- 
ther a power law of energy spectral index q;=1.3±0.2 
or thermal bremsstrahlung with kT=3.9^Q g keV; the 
ridge spectrum appears similar. The features can ac- 
count for most of the hard X-ray flux from IC443 
(§2^). The core of the HXF is marginally resolved 
with the ROSAT HRI, but the low level surround- 
ing emission extends about 10' along the radio-bright 
shell. (§2^, fig. Q). It is spatially coincident with 
the AA20-200 cm fiat spectral {a « 0.2) region (G86; 
fig. 0), and in the general region of, but not uniquely 
coincident with, regions of cloud/shock interactions 
(DeNoyer, 1979; fig. |). It is located outside the 95 
percent confidence error circle of the nearby EGRET 
source (EHKS; fig. ^). No strong periodicity is ob- 
served from the feature, with a pulse fraction upper 
limit of 33% (§p7^), and there is no known pulsar near 
it. 

3.2. Origin of the Emission 

We focus here on the origin of the emission from the 
HXF. Because of the uncertainty of the spectral pa- 
rameters we know substantially less about the ridge. 
While some of our discussion below applies to all the 
emission (e.g. our conclusion that it is non-thermal), 
a more definitive understanding of the nature of the 
ridge emission (which could be different from that of 
the HXF) awaits a better quality X-ray spectrum. 

The data suggest that the most likely mechanism 
for producing hard X-rays is synchrotron emission. 
Before addressing whether this emission arises from a 



plerion or shock accelerated electrons, we first discuss 
our rationale for ruling out alternative mechanisms, 
including bremsstrahlung from a thermal and a non- 
thermal population of electrons, and inverse-Compton 
scattering. 

Thermal Bremsstrahlung A thermal bremsstrah- 
lung model provides an adequate fit to the ASCA and 
the Cinga spectra separately. The inferred tempera- 
tures are discrepant, however (« 15 keV for Cinga 
versus « 4 keV for ASCA). In contrast, a power law 
model yields the same photon index. If the 2-20 keV 
spectrum is better characterized by a power law, then 
the temperature obtained from a thermal model fit 
from each instrument would yield a temperature char- 
acteristic of the instrumental bandpass, which is ex- 
actly the case here. Thus, a thermal bremsstrahlung 
model is inconsistent with the ASCA and CINCA ob- 
servations taken together. 

Inverse Compton One possibility discussed by 
WAHK is that X-rays are produced by inverse-Compton 
scattering of infrared photons off electrons. Gaisser 
et al. (1996) found that an inverse-Compton compo- 
nent, scattering off the microwave background, would 
have a photon index of r=1.5. On the other hand, 
if the scattering photons were locally produced, one 
would expect the hard X-ray emission to arise from 
IR-bright regions. The localization of the hard X- 
ray emission a restricted region, the brightest part 
of which has no IR counterpart argues against an 
inverse-Compton origin. 

Accelerated Bremsstrahlung It has been sug- 
gested that electrons accelerated in the shock to MeV 
energies may be responsible for the generation of 
hard X-rays via bremsstrahlung in supernova rem- 
nants (Sturner et al. 1995; Asvarov et al. 1990). The 
fact that the radio spectral index around the hard X- 
ray feature is flatter than elsewhere in IC 443 is strong 
evidence against this mechanism. If we assume that 
the same acceleration mechanism produces both the 
MeV and the GeV electrons (responsible for the radio 
emission), then there will be fewer MeV electrons in 
the hard region than elsewhere in the remnant. Thus, 
if this mechanism were operating efficiently, this re- 
gion would be dimmer in hard X-rays than every- 
where else in the remnant. 



Synchrotron Radiation The most straightforward 
mechanism for producing a power law X-ray spectrum 
is synchrotron radiation. None of the observations 
contradict this interpretation; all are consistent with 
it. In particular, the fact that the radio spectrum 
of this region is flatter than elsewhere argues that at 
higher energies (including the X-ray band), the syn- 
chrotron flux will be enhanced over elsewhere in the 
SNR. In fact, as we show below, the physical size of 
the HXF corresponds well with the synchrotron loss 
time of X-ray producing electrons. 

WAHK dismissed synchrotron emission as a likely 
model, because an acceleration mechanism efficient 
enough to accelerate electrons to ^10 TeV seemed 
unlikely and they calculated synchrotron loss times to 
be less than the age of the SNR. Their dismissal was 
based on the assumption that the hard X-ray emission 
is spatially uniform; we now know from this ASCA 
observation that the hard emission arises primarily 
from localized regions. 

X-ray synchrotron emission in a SNR can be pro- 
duced by high energy electrons accelerated and in- 
teracting with a magnetic field in one of two loca- 
tions: the magnetosphere of a pulsar (giving rise to 
a plerion) or the forward shock of a SNR. The obser- 
vations support both interpretations to some extent, 
but more issues arise from the presence of a plerion 
than from an isolated region of intense shock acceler- 
ation. 

The compact size of the hard X-ray feature is the 
primary evidence in favor of a plerionic interpretation. 
The lack of X-ray pulsations is not necessarily a prob- 
lem: the upper limit on the pulsed fraction is higher 
than other pulsars like Vela (Pravdo et al. 1976). Nor 
is the existence of extended emission without an obvi- 
ous embedded point source unprecedented: the SNR 
3C 58 has a small extent, a power law spectrum and 
shows no pulsations, but is generally regarded as a 
plerion (Hclfand et al. 1995). 

Difficulties with the interpretation arise when try- 
ing to associate a plerion with IC443. Wilson (1986) 

has shown that the median X-ray:radio ( ,^^,^'11^^ ) 

■' , V 10^-10" Hz/ 
flux ratio for plerions is about 1, with a range from 

about 0.1-500. The flux ratio of the HXF is about 

900, which would give it the highest known X-ray:radio 

flux ratio of any plerion. 

Furthermore, if the feature is a plerion associated 

with the IC 443 shell it requires an extremely high 



projected velocity of 5, 000 (^^) (^) (^j^^ 

kms~^, for distance d, angular distance 9 from the 
explosion center, and age t. IC 443 is in a complex re- 
gion of interstellar space, and the diffuse emission is 
thought to the product of multiple supernova events 
( AA94) . The proper motion problem is circumvented 
if the hard X-ray feature represents the site of the 
most recent explosion which we can take from Chi- 
nese records to have occurred in A.D. 837 (WAHK). 
If that were the case, then it is surprising that this ex- 
plosion has apparently not affected the temperature 
or surface brightness distribution of the diffuse emis- 
sion. The former suggests a more centrally located 
explosion; the latter an explosion in the northeastern 
quadrant of the SNR. 

On the other hand, the presence around the hard 
feature of many interesting structures associated with 
collisions between the shock front and concentrations 
of material suggests the hard X-ray emission is re- 
lated to them. In particular, shock/cloud collisions 
can locally enhance particle acceleration (JK93). We 
consider first whether shock acceleration is a plausible 
source of the hard X-rays, and then how the morphol- 
ogy of the feature might arise. 

Dickel & Milne (1976) found Faraday rotation mea- 
sures in IC 443 of about 200 rad m^^ where the field 
is along the line of sight, and close to zero where it 
is perpendicular. Taking the density and path length 
measurements from the X-ray (PSSW), we estimate 
the magnetic field strength B to be: 



B w 1.23 X 



200radm-2 
Scm^-' X 0.1 pc 



500 mG (1) 



The synchrotron photon/electron energy relation there- 
fore can be written: 



E, 



photon 



5 keVx 



B 



x(: 



Ee 



(2) 



500 AiG '20 TeV 

where ii^photon is the observed photon energy pro- 
duced as synchrotron radiation from an electron at 
energy E^- So while radio emission (i^photon ~ 10~^ 
keV) is produced by GeV electrons, production of 5 
keV X-rays requires the presence of ^^20 TeV elec- 
trons. 

In attempting to explain the non-thermal compo- 
nent dominating the X-ray emission from SN 1006 
as synchrotron emission from highly relativistic elec- 
trons, Reynolds (1996) has shown that it is possible to 
accelerate electrons (and ions) in SNR shocks to ener- 
gies exceeding 100 TeV. While his model is not strictly 



applicable here, some aspects of it can be used to es- 
tablish the plausibility of shock accelerated electrons 
as the source of the hard X-ray emission in IC443. 

In the simplest models, the synchrotron spectrum 
is expected to be a broken power law. The vari- 
ous break frequencies correspond to electron energies 
where two canonical times scales equate, such as the 
synchrotron loss time and the acceleration time scale. 
At each break, the spectral index increases by about 
0.5. The electron diffusion time is also important in 
determining the steepening of the spectrum, but its 
influence is less straightforward to estimate. 

In the IC 443 hard feature, the difference between 
the radio and the X-ray spectral index is approxi- 
mately 1.1, suggesting the need for at least two spec- 
tral breaks between the bands. The time it takes for 
an electron of energy E to lose half its energy via syn- 
chrotron radiation in a magnetic field of strength B is 
given by: 



/SOO/iG 



20 TeV\ 



X 2.5 years , (3) 



which is approximately the light crossing time of the 
HXF. Using Equations 13 and 0, and assuming a SNR 
age between 1000 and 5000 yr (WAHK, AA94), we 
find that the break frequency where the synchrotron 
loss time becomes comparable to the age of the SNR 
occurs in the far infrared. Additionally, by assuming 
a shock velocity on the order of 1000 km s""'^, we find 
the break frequency associated with the equivalence 
between the acceleration and loss times is in the hard 
ultraviolet band. The frequencies at which these two 
spectral breaks occur supports the idea that the flat 
radio and the hard X-ray spectra are both produced 
by shock acceleration. 

While Reynolds' (1996) formalism facilitates the 
plausibility argument above and is especially sugges- 
tive with regard to the ridge region, his models for 
producing X-rays from shock accelerated electrons 
from a large segment of the SNR shell clearly do not 
apply to the HXF in IC443. A different mechanism 
must be operating there. Such a mechanism has been 
suggested by JK93. 

Figure || shows BGBW's H2 image of shocked 
molecular gas. This shows that the molecular cloud 
has a clumpy ring structure exterior to the site of the 
supernova event. Dickman et al. (1992) found that if 
the ring is circular, it must be inclined ~ 51° from the 
line of sight and have a diameter of ~ 9 pc. Doppler 
shift information shows the southern edge of the ring 



(near our HXF) expanding toward us, and the north- 
ern edge away ~ implying that the ring is now expand- 
ing at a deprojected velocity of 25 km/s. While, we 
agree with this picture in principle, we point out that 
to second order the ring is not circular; most notice- 
ably there is a section of negative curvature surround- 
ing the HXF (Fig. 0). It is this change of curvature 
that may create the HXF. 

JK93 simulated cosmic-ray acceleration as shocks 
impact dense clouds - which is the case here. They 
found that after a shock/cloud impact, particle accel- 
eration is significantly enhanced in the "tail shock" 
emanating from the cloud in the direction of shock 
propagation. Their results suggest that each clump in 
the molecular ring could produce enhanced shock ac- 
celeration in a long tail behind it. Because of the neg- 
ative curvature of one section of the ring, the flow will 
be refracted, thereby "focusing" the enhanced shock 
acceleration at the hard X-ray feature where the tail 
shocks overlap. This interaction between these tail 
shocks will both increase the effective shock veloci- 
ties and provide a complicated structure which will 
require further study. 

The total hard flux from the diffuse emission along 
the shell (the ridge) is about half that from the HXF, 
but the ridge may extend out of the field of view. It 
corresponds to an arc bright in H2, and radio con- 
tinuum. While we cannot yet be certain that the 
emission in non-thermal, the higher ambient density 
as indicated by the infrared line emission argues that 
the temperature of shock-heated material should be 
lower than elsewhere in IC 443 - not higher. The ar- 
guments made above for the HXF against emission 
mechanisms other than synchrotron hold for the ridge 
region as well. While we withhold judgment until a 
higher quality spectrum can be obtained, we suggest 
that we are seeing evidence here as well for TeV elec- 
trons. Whether these arise from the JK93 mechanism 
or are more closely related to the forward shock is still 
an open question. 

4. Conclusion 

In this paper we have presented evidence for a lo- 
calized region of particle acceleration to electron en- 
ergies over 20 TeV, within the shell type supernova 
remnant IC 443. Because of the unique positions of 
the hard X-ray feature and shocked clouds, we be- 
lieve that the HXF is caused by the focusing of tail 
shocks around each molecular cloud, thus enhancing 



particle acceleration in this location. In addition we 
have found a ridge of hard emission which in co- 
incident with both the radio-synchrotron shell and 
shocked molecular gas, but whose spectral parame- 
ters are poorly constrained. 

More observational and theoretical research is needed 
to fully understand this phenomenon. Multi-frequency 
and polarization observations with the VLA will allow 
us to both measure the feature's spectral index and 
look for other indicators of on-going particle accel- 
eration. A deep single-dish radio observation of this 
area is necessary to verify that the HXF is indeed 
non-plerionic. The most important future observa- 
tion will be to observe the full extent of the ridge and 
the shocked molecular gas with ASCA to better char- 
acterize the spectrum and image the full extent of the 
hard non-thermal emission. 

On the theoretical side, this discovery provides 
a nice opportunity to model particle acceleration 
and cloud/shock interactions with significant obser- 
vational constraints. 

We thank the following people for contributing to 
these results: Richard Mushotzky for bringing the 
hard feature to our attention; Michael Burton for 
kindly providing his H2 emission line image in dig- 
itized form; Steve Reynolds and Tom Jones for their 
insights into shock acceleration; Matthew Baring and 
Ocker de Jager for their thoughtful comments and 
suggestions. 
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